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Abstract: In this work, we study the large scale structure formation in the modified 
gravity in the framework of Palatini formalism and compare the results with the equivalent 
smooth dark energy models as a tool to distinguish between these models. Through the 
inverse method, we reconstruct the dynamics of universe, modified gravity action and 
the structure formation indicators like the screened mass function and gravitational slip 
parameter. Consequently, we extract the matter density power spectrum for these two 
models in the linear regime and show that the modified gravity and dark energy models 
predictions are slightly different from each other at large scales. It is also shown that the 
growth index in the modified gravity unlike to the dark energy models is a scale dependent 
parameter. We also compare the results with those from the modified gravity in the metric 
formalism. The modification on the structure formation can also change the CMB spectrum 
at large scales however due to the cosmic variance it is hard to detect this signature. We 
show that a large number of SNIa data in the order of 2000 will enable us to reconstruct 
the modified gravity action with a suitable confidence level and test the cosmic acceleration 
models by the structure formation. 
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1. Introduction 

For more than a decade, the positive acceleration of the universe is one of the challenging 
questions in physics ||l|, |2|. The physics and the mechanism behind this phenomenon is 
unknown. The Cosmological Constant (CC) is the most straightforward suggestion to 
describe an accelerating universe However the fundamental questions like the fine 

tuning and coincidence problems opened new horizons to introduce the alternative models 
g like Dark Energy (DE) g and models of the Modified Gravity (MG) 0. The other 
motivation for the modified gravity models is unifying the dark matter and the dark energy 
problems in a unified formalism ||8|. 

The MG theories need to be tested in two domains of the observations: (a) cosmological 
scales 1^] and (b) local gravity [10|. The simplest modification of gravity is the extension 



of Einstein-Hilbert action where instead of Ricci scalar, we use a generic function of f{R)- 
These modified gravity theories are divided into the Metric and Palatini formalism |Tl| . 
There are a lot of debates in the literature on the viability of this type of modified gravity 



model concerning the local gravity tests |12] where a vast range of them had been falsified 



|13, Also many models are proposed to evade the falsification criteria |15, T^, 17]. The 
Palatini formalism is a second order differential equation which has more simpler structure 
than the fourth order metric formalism. Since all the conventional field equations are 
described by the second order differential equation, it seems that gravity should also follow 
this rule. So this is the main advantage of working with the Palatini formalism. From 
cosmological point of view, we will show that for a given dynamics of the Universe there 
is a one-to-one map to the corresponding action in the Palatini formalism. Also in this 
formalism we don't have the instability problem as in the metric formalism. While there 
are advantages in this formalism, there is challenging questions in the Palatini formalism 
as the microscopic behavior of the matter fields where in the Einstein frame, these models 
disagrees with the standard scalar-matter coupling theories Means that we can not 
apply the perturbation theory in these scales as the gravity field is directly relates to the 
energy-momentum tensor, in contrast to the Einstein gravity where the amount of the 
perturbation of the metric at a given point is averaged over all the space. Although there 
are debates on how we should do averaging procedure over the microscopic scales [^], it 
has been shown that the problematic microscopic behavior of this theory can be solved in 
the case of f{R) very near to the ACDM model [11|. 

Keeping in mind that there are many fundamental questions in the Palatini formalism 
, here in this work we mainly focus on the behaviors of the modified gravity models at the 
cosmological scales and their predictions on the structure formation. The main question we 
address in this work is the investigation of cosmological probes in the structure formation 
to distinguish between CC, DE and MG models |20|. We first reconstruct the dynamics of 
the background and the corresponding MG action via the inverse method from the SNIa 
data, described in |21]. The dynamics of the background (i.e Hubble parameter H = H{z)) 
can only distinguish between CC and the alternative models |21| however in the level of 
probing the expansion history of the Universe, MG and DE models are not distinguishable. 
Consequently, we use the structure formation probes as the promising tool to distinguish 



between alternative models ||2^, 23]. It is worth to mention that the alternative models can 



be extended to a complicated forms such as interacting dark energy-dark matter models 
]24] or clustering DE theories, which can not be distinguished from the MG models even 
in the level of the structure formation probes [p5|]. 

In this work we focus on the structure formation issue in the modified gravity-Palatini 
formalism and compare the observational effects with that in the smooth dark energy 
models (sDE). We derive the screened mass function defined as the fraction of effective 
gravitational constant to the standard gravitational theory as, Q = G(,f ^ jG^ . The effective 
gravitational constant appears in the Poisson equation depends on the scale of the structure. 
The other relevant parameter appear in our calculation is the gravitational slip parameter 
7 = — ^>/^', which is the fraction of spatial perturbation of the metric to the perturbation 
of the time-time element. The screened mass function as well as the gravitational slip 
parameter in DE models are equal to one in GR while in MG models they are time and scale 
dependent parameters . This scale dependance of structure formation in modified gravity 
theories is a well known effect discussed in literature ]16, 17, 26]. A combination of screened 
mass and gravitational slip parameters recently have been used for distinguishing the MG 



in metric formalism with the alternative models via analyzing the structure formation and 
weak lensing |27]. We show that the derived reconstructed matter density power spectrum 
and growth index in this model are affected by the screened mass function. Also we study 
the effect of gravitational slip parameter on the Integrated Sachs- Wolf (ISW) and compare 
it with recent observations of WMAP data. As a comparison with the metric formalism, 
we compare our results with that of structure formation in the metric formalism. 

The structure of this article is as follows: In section ^ we review the standard equa- 
tions governing the relativistic structure formation. In section ^ we re-derive the structure 
formation equations for DE and Palatini MG models. Also in this section the metric for- 
malism structure formation equations are derived for comparison with that in the Palatini 
formalism. In section Q we extract the screened mass function and gravitational slip pa- 
rameters, reconstructed from a modified gravity equivalent to a dark energy model and 
the effect of screened mass function on the power spectrum of the structures as well as 
the growth index is studied. Also the screened mass in the metric formalism is obtained 
in order to show the scale dependence of structure's growth. Details of the calculation in 
metric formalism is given in Appendix A. In section |5| we investigate the effect of gravi- 
tational slip parameter on CMB power spectrum due to the ISW effect and compare the 
results with recent WMAP data. The conclusion is presented in section ^ 



2. Structure Formation: Perturbation theory 

The relativistic structure formation in the universe can be studied by the perturbation 
of the FRW metric and homogenous energy-momentum tensor of the cosmic fluid . A 
standard technic is decomposition of the metric and energy momentum tensor perturbations 
into the scalar, vector and tensor modes where each mode evolves independently |2^]. Since 
in this work we only deal with the density contrast evolution of the structures, we study 
the evolution of the scalar modes in our study. One of the ambiguities in the perturbation 
of metric is the gauge freedom, similar to what appears in the theory of Electromagnetism. 
In fact the problem arises due to the excess number of the degrees of freedom (i.e. number 
of variables) compare to the number of independent field equations. We fix this freedom 
with so-called Newtonian gauge as follows [30| ]: 

ds'^ = -[1 + 2^{x,t)]dt'^ + a\t)[l + 2<!>{x,t)]glf . (2.1) 

It should be noted that the two scalar perturbations of ^ and ^ are gauge invariant 
variables and hence observable parameters, as the electric and magnetic fields in the elec- 
tromagnetism [^]. In order to interpret the physical meaning of the perturbations, let us 
look to the the Newtonian limit. In the Newtonian limit, the geodesic equation reduces to 
X = — V^'(x, t). On the other hand the Einstein equation in the weak field regime reduces 
to the Poisson equation as follows = AnGpS. 

The dynamics of the perturbations represent the evolution of the structures which can 
be obtained from the Boltzman and the Einstein field equations. First we consider the 
equations arise from the perturbed field equations. The (0-0) component of perturbed 



Einstein equation {i.e.6Gu = uSTy) in the fourier mode is given by |29| 



k^^ + 3H{^ - H'^) = 4TTGa'^p6, (2.2) 

where derivative is in terms of the coordinate time defined in the metric and we ignore 
the subscript k for simphcity. On the other hand the spatial component of the Einstein 
equation, after applying the projection operator of kik^ — l/Sdj in Fourier space, results 
in: 

A:2($ + $) = -327rG;0°e2, (2.3) 

where the right hand side of this equation is derived from the first order perturbation of 
the distribution function of the cosmic radiation. Here ©2 is the the second moment of 
anisotropic stress of radiation where l-st moment is defined by 

(T^ y"^'^P/(M)e(M), (2.4) 

and Pi{n) is the Legendre polynomials. Q{fJ-) is the temperature contrast of radiation and 
defined by: 

T{x,f,t) = To[l + e{x,f,t)], (2.5) 

where Tq is the mean temperature of radiation, and the fj, is the cosine of the angle between 
the momentum of a fluid and the direction that observer looks at the temperature contrast. 

In the matter dominated universe where the energy density of radiation is negligible, 
the right hand side of equation ( p.3p is zero, hence ^ = — In the next section we 
will show that for the modified gravity models even in the matter dominant epoch the 
scalar perturbations have no more simple relation as in the general relativity. In order 
to discriminate the modified gravity models from sDE, we use the new parameter of 7, 
so-called gravitational slip parameter as [^ : 

7 = -|, (2.6) 

In the next section we relate the slip parameter to the observable parameters in the struc- 
ture formation. Any deviation from 7=1 can show the crackdown in our assumptions. 

As we noted before, the other set of the equations governing the structure formation 
are the perturbed Boltzman equations. For a given perfect fluid for the universe with the 
energy momentum tensor of T^jy = {p + p)u^Uu + POfiu-, we can obtain the energy and 
momentum conservation equations by T^'^-u = for ^ = and fj, = i as below: 

3<^-3H{^-io)——, (2.7) 



1 + u! a 6p 1 + CO 

e = -„ii-^)e--^e^'Jf-^^^*^, (2.8) 

where 9 is the divergence of peculiar velocity 9 = V ■ v and u is the equation of state. 
Since we assume that the cosmic fluids are already decoupled from each other at the 



late time universe, equations ( |2.7] ) and ( |2.8| ) can be written for each component with the 
corresponding equation of states. 

The goal of structure formation theory is to find the dynamics of the four perturbed 
quantities of ^> and ^ from the metric and 6 and 9 from the energy momentum tensor 
through the field and Boltzman Equations 



3. Structure formation in accelerating universe 

In this section we reexpress the perturbation equations for the structure formation in an 
accelerating universe. We note that the derivation for sub-horizon scalar perturbations in 



the general case of scalar-tensor gravities were obtained in [34| where the Palatini formalism 
is a special case of it. In the first subsection we will introduce the equation of the structure 
formation for a universe filled with a smooth dark energy with an equation of state lode = 
< — ^. In the second part we derive the equation of the structure formation in the 
Palatini MG model equivalent to a Dark Energy model. 

3.1 Smooth Dark Energy Models 

In this subsection, we study the structure formation assuming a DE fills the universe. In 
a generic case, if we assume a clustering dark energy model, the perturbation of Einstein 
equation in the non-radiation dominate epoch can be written as [ p5[ : 



a 



+ 2,H{^ - H^) = iTTG[p^5^ + 5dePde] (3.1) 



A;2($ + ^r) = -327rGa, (3.2) 
where a is the anisotropic stress of DE model. This quantity is a generic property of a 



fluid with flnite shear viscous coefficients [36|. In case of our study, in which we consider 



a perfect fluid for the matter and DE, the anisotropic stress is vanished (i.e a = 0) and 
$ = Using the trace of perturbed Einstein equation 6R = —87rG{36p — 5p) and 



combining with Eq. (3.1), we can eliminate the density perturbation and write a second 



order differential equation for $ as follows p7| ]: 

^ + m^ + {2H + ^H'^)^ = -A-KG5pDE- (3.3) 

In the simple case where we don't have clustering of the DE we set 5pDE ~ and the 
dynamics of $ can be obtained from (p.3|). We can substitute the numerical value of <I> in 



(3.1) and obtain the evolution of the density contrast for the matter as: 

2 



1.2 

— $ + 3/?^$ + 3F$ = 47rG5/9m. (3.4) 



It is worth to mention that DE models with the clustering properties can modify the 



structure formation and consequently change the matter power spectrum |25]. 

In order to compare our results in the following sections with the observation, we 
reexpress the differential equation (^^) in terms of the redshift: 



where E{z) = H/Hq is a dimensionless Hubble parameter, Hq is the Hubble parameter at 
the present time and k = k/ Hq is the dimensionless wavenumber of the structures and JlJ^ 
is the matter density parameter at the present time. 

3.2 f{R) Modified gravity theories in Palatini Formalism 

We take /(i?)-modified gravity in the Palatini formalism as an alternative model for the 
acceleration of the universe. In this formalism the connection and metric are independent 
fields. Variation of action with respect to these fields, results in a set of second order 
differential equation for the metric p8| . 

The action of MG in palatini formalism is given by: 

S[f;g,r,^rn] = -^l d'^xV^fiR) + S^[g^„^^], (3.6) 

where k = SttG and Smig^^ui ^m] is the matter action depends on the metric Qf^^ and the 
matter fields R = 9^'^ R^iuO^) is the generalized Ricci scalar and R^i, is the Ricci tensor 
made of affine connection. Varying the action with respect to the metric results in 

f'{R)R^,{r) - = f^T^u, (3.7) 

where prime is the derivative with respect to the Ricci scalar and T^,/ is the energy mo- 
mentum tensor defined by 

The trace of the field equation results in: 

RF{R) - 2f{R) = kT, , (3.9) 

where T = g^'^T^^ = —p + 3p, and F{R) = df /dR. First of all we want to know the 
background evolution of the Universe in Palatini formalism. We use FRW metric in flat 
universe (namely A; = 0). 

ds'^ = -dt^ + a{tf5ijdx'dx\ (3.10) 

and assume that universe is filled with a perfect fluid with the energy-momentum tensor 
of T'^ ^ = diag{—p,p,p,p), the generalized FRW equations is as below 



r2 _ 1 3(1 + oj)f - (1 + 3u)Rf' 
6(1 - 3a;)/' 



1 + 2 (1 + W) 

In the matter dominated epoch where w = 0, the Hubble parameter reduces to: 



1 , 3 f"(2f-Rf') 
i -r 2 f(Rf"-f') 



In order to study the structure formation in MG models, we perturb the metric in 
the same way as done in the previous part. In the modified Einstein equation, we move 



the extra geometrical terms out of the Einstein tensor to the right hand side of equation, 
recalhng it by j''^^'^"'''^)^ ^he perturbed modified Einstein equation can be written as 
fohows: 



The (0,0) component of the field equation is |4l[] : 



F ■ 2k^ 1 SF^ 

3(2i? + -){H^ _ _ _ci> + _(_ + 3HF)^ 



1 

F 



{-k6p + 5F[3H'^ 



3F' 
4F2 



R k\ 



+ (^l + 3H)5F}, 



(3.13) 



(3.14) 



where 5F is the variation of the F{R) due to the perturbation of the metric. Eq.(3.14) is 
a modified Poisson equation where for F = 1 we recover equation (3.1). G^-component of 
field equation also is given as below: 



1 

2F 



6F-{H + ^)6F -F^ + KpJ^ 



(3.15) 



where = k'-T^i/ip + p) is a gauge invariant quantity. Now we combine equation ( |3.14| ) 
and ( |3.15| ) to obtained the Poisson equation in the modified gravity as |42]: 



k'^(^ = 47rGa^AmpQ{k,z), 



(3.16) 



where 



A^ = 6n + 3H9k/k^ 



(3.17) 



is a gauge invariant quantity in terms of the density contrast |43] and Q{k,z) is so-called 
screened mass. This term in the Palatini formalism is given as 



1 mR-^D, 

+ n ' 

r 1 — m 



(3.18) 



where m = R-^ In F{R) is a dimensionless deviation parameter and D is defined as below: 



D = 6H' + -{ 



3.F. 



2'F' 



+ 6H 



.F R 



F 



2, a 



In the modified gravity theories where the action is close to the Einstein Hilbert (i.e m <C 1), 
the screened mass reduces to [40|: 



Q{k,z) 



F^^ l-m' 



(3.19) 



where ^ = mk'^a~'^R~^ . 

We note that in the definition of the gauge invariant density contrast in Eq.( 3.17| ), 
the extra term added to 5 is in the order of VHVp/(?-, where vh is the velocity of Hubble 



expansion at the size of the structure and Vp is the corresponding pecuUar velocity. In 
order to have an estimation from this term, we take Harisson-Zeldovich spectrum for the 
density contrast. Then the density contrast of a structure changes as 5k = S enter k'^ / H'^ , 
where 6enter = 10~^ is the density contrast of the structure with the horizon size at the 
present time. On the other hand the pecuhar velocity in the linear regime is in the order of 
Vk = Hodk/k. Substituting 5k from the Harisson-Zeldovich spectrum, the peculiar velocity 
obtained as Vk = k/HQ5enter, or VkVn — 10~^. For sub- horizon scales where 5^ ^ 10~^, 



we can neglect the second term in equation ( 3.17 ). The results within this approximation 
is identical to the case where we choose comoving gauge in the analysis. In this gauge, the 
peculiar velocity is set to zero [|^. 

Another indication of the modified gravity is that, the right hand side of the Eq. ( p.3| ) 
in the late time universe is non-zero. Using the projection operator of k^k^ — 1/3(5| on the 
perturbed field equation of ( 3.13| ) results in 



4C, 44 



= $ + ^r = 



5£ 
¥ 



(3.20) 



where 5F represents the perturbation of the action due to the perturbation of the metric 
which can be rewritten as 5F = Fji5R = dF/dR 5R. We note that for any deviation of the 
gravity law from the Einstein-Hilbert action in Eq. ( 3.20| ), <I> + ^' 7^ and the screened mass 
function is not equal to one. This may affect on the growth of the large scale structures 
and consequently modify the expected power spectrum of the structures. 

The deviation of the Ricci scalar 5R due to the perturbation of the energy-momentum 
tensor is given by the trace of field equation (^^) as below: 

8ttG5T 



5R 



RF 



(3.21) 



Where 5T is trace of perturbed energy-momentum tensor and in the case of matter domi- 
nated era is given by 5T = —5p. We substitute ( 3.21 ) in ( 3.20 ) and change the derivatives 
from the Ricci scalar to redshift, so the Ck in Fourier mode is given by: 



a = ^ + * 



F] 



(3.22) 



where Jim is the matter density parameter in redshift oi z, R = R/Hq and 5rn{k, z) is the 
matter density contrast which depends on wavenumber and redshift. 

Now we can define the gravitational slip parameter in terms of MG free parameter F. 
Using Eqs. ( p.6| ) and ( |3.20| ), the slip parameter can be written as: 

Fr5R 



l{k,z) 



1 + 



^F 



(3.23) 



On the other hand we want to write the gravitational slip parameter independent of the 
matter density contrast. The combination of Eqs. (|3.16| ), (|3.21 ) and ( 3.23| ) results in: 



i{Kz) = [i 



FQ(A:,z)(i?f (f)-i 



FY 



(3.24) 



We can also rewrite the gravitational slip parameter in terms of deviation parameter as 
below: 

2m{l + zfk"^ 



7 = [1 



RFQ{m - 1) 



(3.25) 



In the following sections we will apply the gravitational slip parameter in the Integrated 
Sachs- Wolfe effect (ISW) and study the deviation of the CMB power spectrum from the 
dark energy models. 

In order to find the evolution of the matter density contrast in MG, We start with 
the conservation of energy momentum tensor in Eqs. ( |2.7| ) and ( |2.§| ). For the matter 
component of the perturbation a; = and for the scales larger than Jeans length 6p = 0. 
The energy-momentum conservation equations reduce to: 



1 

--( 

a 



3«>, 



(3.26) 



-He + —^f. 

a 



(3.27) 



Now in order to find the evolution of the density contrast, we take the time derivative from 
Eq.( |3.26| ) and eliminate 9 by substituting it from equation (3.27). The density contrast is 
obtained in terms of the potentials ^' , <I> as below: 



■ 

5 + 2H6 + -^^ 



-6H^ - 3$ 



(3.28) 



In order to determine the evolution of density contrast, we apply the modified Poisson 
equation, Eq.( 3.1(^ ) and the gravitational slip parameter, Eq.(2^) in Eq.( p.28 ). We find 
the differential equation as below: 



[1 +3^Q]Sm + [2H + 6^Q]5r, 



(3.29) 



S^iH^Q + HQ - Q)]6„ 



0, 



where T = AttGp^. Eq.( p.29| ) is the general form of the differential equation for the 
evolution of matter density. In the limit of quasi-static approximation <I> ~ with the 
assumption of studying the scales deep inside the horizon, where ^ 



'H\2 
^ k > 



0, Eq.(|33g|) 



reduces to the equation defined as [43|: 



5^ + 2H5m - A-KGpmSml ^Q{k, z) = 0. 



(3.30) 



For the case of GR, we will have ^ = Q = \ and we recover the standard equation for the 
evolution of density contrast [40|. 

In the following sections we will solve this differential equation for a suitable modified 
gravity, reconstructed from the SNIa data. 



3.3 f{R) Modified gravity theories in Metric Formalism 

As we mentioned earlier there are two different approaches for the modified gravity models. 
It seems that one of the possible tools to distinguish them is the evolution of the large 
scale structures in these two formalism. Here in this section we study the the growth of 
structures in the metric formalism and obtain the evolution of density contrast and the 
screened mass function. Further investigation in metric formalism for the non-linear regime 
of the structure formation will be done in our future work. 

The action of the modified gravity in metric formalism is given by: 



S = j <fx^f{R) + SmiOf^u, (3.31) 



where the connections V^^ are usual connections defined in terms of the metric g^^ . The 
field equations obtained by varying the action (3.31) with respect to g^^: 



F{R)R^,{g) - ^f{Rg^,u - V^V,F(i?) + gf^^nFiR) = SttGT^, (3.32) 
and its trace is 

3nF(i?) + F{R)R - 2/ = SvrGT. (3.33) 

The evolution of density contrast will be derived from perturbed Einstein equations in 
deep inside the Hubble radius is (for more details see Appendix A): 



3H5F - (sH + 3H^ - 6F - 3HF^ - 3F{H'^ - 6) - K6p) 

(3.34) 



6F + 3H5'F + (^ - ^)SF = '^6pm + F{3H^ + ^ - 36) + (2F + 3HF)^ - \f5R, (3.35) 
3 3 3 



5F 

m + ^ = — -. (3.36) 
F 

In the quasi-static regime where \F\ <C \HF\ and \F\ < \H'^F\, the dominate terms will be 
k'^/a'^ and Kp and with the assumption of small deviation from GR the equation governing 
the evolution of matter density will reduce to [ 43 1 : 



Sm + 2H6m - AnGj^l^QmetSpm ~ 0, (3.37) 

where jmet is the gravitational slip parameter and Qmet is the screened mass function, both 
defined in the metric formalism. The screened mass function in the metric formalism is 
given by: 



Q = Geff/G 



4 1 M^a 



2^2 



(3.38) 



3 3A;2 + M2a2 

where = ^{^)~^, which in the case of Af — t- oo we will recover the GR case Q = 1. 
The equation above indicates that in large scales k <C Ma, we will recover Qmet — ^ 1; while 
in the small scales k ^ Ma, the screened mass function converge to 4/3. Meanwhile by 
assumption of the same background dynamics, in the Palatini formalism the screen mass 
function from Eq.( 3.18| ) predicts a different scale dependance. Consequently the structure 



formation observables will have different values for metric and Palatini formalisms. It is 
crucial point to indicate that for this comparison, we should also extract the correspond- 
ing modified gravity action in the metric formalism. For this task we refer to modified 
Friedmann Equations in metric formalism obtained from field equation ( |3.32| ): 

3FH^ = SnGp^ + ~ ^ - 3HF (3.39) 
-2FH = 87rGpm + F - HF 

where dot represent derivative with respect to physical time which can be reexpressed in 
terms of redshift as: 



F + 



2 E 
1 + z "B 



2F E' 31^^(1 + z)_ ^3^^Q^ 



l + z E ^2 



By substituting the reconstructed Hubble parameter E' = in Eq ( |3.40 ) from the equiv- 



alent sDE model and solving the differential equation with GR boundary conditions in 
higher redshifts (i.e. z ~ 100), we can find F = F{z) and eliminating the Ricci scalar 
derive the action [^] . We will review the procedure of deriving the equivalent action with 
a given dark energy model in Palatini formalism in the following section. 

4. Structure Formation: Modified gravity versus Dark energy 

In this section we first review the dynamical equivalence of a modified gravity with a smooth 
dark energy model. Knowing the dynamics of the universe from the observational data, 
we extract the appropriate action for the modified gravity as described in a recent work by 
Baghram and Rahvar From the action, we obtain the screened mass function and the 
gravitational slip parameter. Then we calculate the evolution of the large scale structures 
for the MG and compare it with that from the dark energy equivalence, assuming the same 
initial condition for the density contrast. Finally we compare the power spectrum of the 
structures as a discriminating tool in these two models. 

4.1 Geometrical Equivalence of the Modified gravity theories with the smooth 
Dark Energy models 

In this part, we review the equivalence of dark energy and modified gravity in the back- 
ground dynamics of the universe (i.e. H = H{z)) In the standard GR, from the 
Friedman equations the energy density is related to the dynamics of the background as: 

= f P, (4.1) 

a K, 
a b 

where p and p are the overall energy density and the pressure of the cosmic fiuid. The 
effective equation of the state for the cosmic fiuid can be written as: 



(4.2) 



Separating the contribution of the dark matter and dark energy in the matter budget of 
the universe, the equation of state of the dark energy is related to the effective equation of 
state as follows: 



E^{z)uJeff 



(4.3) 



where is the matter density parameter at the present time. For a given modified gravity 



model, we can substitute the Hubble parameter in equations (|4.2| ) and ( [4^ ) to obtain the 
equivalent equation of state of a dark energy model. It should be noted that the geometrical 
equivalence between the two models is a one-to-one relation. 




Figure 1: The dimensionless Hubble parameter, E{z) = H{z)/Hq versus redshift derived from 



solving the Eq.(4.3) for a dark energy model with the equation of state defined in (4.7). 



On the other hand, for a given equation of state of a dark energy model, ujde{z), one 
can extract the dynamics of the universe, E{z). Knowing the dynamics of the universe, 
from Eq. ( 3.12| ) , the Hubble parameter is given in terms of the derivatives of the action and 
the Ricci scalar. The Ricci scalar in Palatini formalism is given by: 

3V„V"F 



F 



in which R{g) is the Ricci scalar in terms of the metric and R{g) = 6H + 12H^ . 
FRW metric, the Ricci scalar in the Palatini formalism is given by |p5[| : 

F 



R 



-QHH'{1 + z) + 12H^ - -H^{1 + zfi^f 



+ 3(1 + zfHH'^ - 6H^{1 + z)^ 
F" 

+ ZH'^{l + zf — . 
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For the 
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Figure 2: The deviation parameter <^ = df /dR — 1 versus redshift derived from solving the differ- 
ential equation (4.6). Here the modified gravity is obtained from the dynamical equivalance of a 
dark energy with an equation of state given in (4.7). 



Now we substitute the Ricci scalar in Eq.( |3.12j ) and with some algebraic manipulation we 
obtain a differential equation for the evolution of F in terms of redshift: 

F" — + F'{— + 



2 F 1 + 

where prime represents the derivative in terms of the redshift. Assuming p = and p = 
po{l + z)^, we eliminate the density in favor of the Hubble constant, using npo = 3HqQ^. 
Finally we can solve the differential equation for a given Hubble parameter and obtain F 
in terms of redshift. 

Here we choose a smooth dark energy ansatz with the equation of state of [^]: 

UJ = Uo + UJa— — , (4.7) 
1 + z 

where ujq and coa are the free parameters of the model. In the case of = —1 and = 
we recover the ACDM model. Here in our simulation we choose loq = —1.2, oja = 1.5 and 



use Eqs. (4^) and ( [4. 3D to extract the corresponding Hubble parameter for this equation 
of state as shown in Fig.(|l|). On the other hand substituting the Hubble parameter in Eq. 
( |4.6| ) results in F in terms of redshift. We set the boundary conditions in redshift of z ~ 100, 
where we anticipate to recover the general relativity condition (i.e. F{z = 100) = 1 and 
F {z = 100) = 0). Applying the Hubble parameter in Eq.([4.5D we get Ricci scalar in terms 



of the redshift. On the other hand we substitute the Hubble parameter in (4.6) and get 




Figure 3: The modified gravity action f{R), which is obtained from integration of deviation 
parameter ^ 

F in terms of the redshift. We plot the parameter of the action deviation, ^ = F — 1 as 
a function of z in Fig. (^). Ehminating the redshift between the F and R results in the 
derivative of the action in terms of the Ricci scalar, F = F(R). One step further we do 
integration and find the modified gravity action f{R) in terms of R, as shown in Fig.(^). 
A numerical fit to this function is: 

f{R) R-ae'^, (4.8) 

where a ~ 0.1 and Rq ~ 8.48i^Q is the Ricci scalar at the present time. It is obvious from 
the form of Eg. ([4.^) that for large scalar Ricci we will have the Einstein-Hilbert action. 

It is worth to mention that for a smooth potentials of the scalar fields in the quintessence 
models, Bassett et al. (2002) and Corasaniti et al. (2003) [^] show that the equivalence 
dark energy model with this class of the quintessence models can be identify with two free 
parameters of wq and the derivation of the w at the present time (i.e. w{z) = wq + f+f )• 
We did our analysis based on this class of the dark energy models and didn't test in the 
cases that w{z) has dramatic variation with to redshift. 

4.2 Reconstruction of the Action in the Modified Gravity 

Our approach for reconstruction of the action is based on the measurement of the distance 
modulus of the SNIa data. From the SNIa data we can make a continues function for 
the distance modulus as a function of redshift (i.e. fJ-{z)), using the smoothing method 
1 48, 49, 50] as described in pi]]. The smoothing method suffers from the low statistics of 
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Figure 4: The reconstructed Hubble parameter with la confidence level obtained from the recon- 
struction of 100 ensemble of SNAP data. 
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Figure 5: The reconstructed deviation parameter ^ = df /dR— 1 with la confidence level obtained 
from the synthetic SNAP data. The uncertainly in ^ indicated with the solid lines, represents 
different realizations of SNIa data in an ensemble of 100 data set. 



the SNIa data where it has been shown in that at least ~ 1500 SNIa data with the 
quality of Union sample data set |52| is essential to generate a reliable distance modulus. 

Here we repeat the same procedure with 1998 synthetic SNIa data with the quality 
predicted by SNAP project in the redshift range of 0.1 < z < 1.7. We also add 25 SNIa 



data from Nearby Supernova factory project in the redshift of z < 0.1 to have a full 
range coverage of SNIa data Now we use the hypothetical cosmological model as 

described in Eq. ( [4.7D for generating the corresponding luminosity distance for the SNIa 
data with a systematic error of Csys = 0.07 in SNAP. 

After generating the synthetic distance modulus, we obtain the best curve resulting 
from the smoothing of the simulated data. For the smoothing discrete data, we first take 
an initial guess model for the distance modulus based on AC DM model. Then we smooth 
the subtraction of the distance modulus of simulated data from the guess model using a 
Gaussian function as below: 



Az) = N{z) Y,[i.^{z) - exp ^\ , (4.9) 



where 



iV(z)-i = J^expi^^^], (4.10) 



Zi represents the redshift of each SNIa in the simulated SNAP sample. The sum term is 
considered for all 2023 SNAP simulated data sample, /i'^*™(zj) is the simulated distance 
modulus and fJ-^{z) is a continues function generated by the guess model. N{z) is a normal- 
ization factor, A is a suitable redshift window function and n^{z) is a smoothed continues 
function for the residual of luminosity distance in terms of redshift. Finally we add the 
smoothed distance modulus to the guess model to generate the new smooth function for 
the distance modulus: 

M^(z)=^^(z)+/(z). (4.11) 

We repeat this procedure using fj,'^{z) as the new guess model. It can be shown that after 
a finite time of irritations, the of the smoothed function with respect to observed data 



will converge to a fixed value[51|. Now we have a continues distance modulus function in 



which the Hubble parameter in terms of redshift can be obtained by: 

H{z) = iU^r'. (4.12) 
az 1 + z 

In order to calculate the uncertainty of this procedure, we generate an ensemble of the SNIa 
data and extract the corresponding Hubble parameters. Fig.(^ shows the reconstructed 
Hubble parameter with la level of confidence. 

Now we can extract the modified gravity action in Palatini formalism through the 



method of inverse problem |53] in the Palatini formalism as discussed in the last part 



Applying the Hubble parameter from Eq.( 4.12| ) in Eq.( [4.6| ), we obtain ^ in terms of redshift 



as plotted in Fig. (^) with la level of confidence. Numerical integration of F{R) provides 
the action, f{R). 

4.3 Reconstruction of structure formation indicators and matter density Povi^er 
Spectrum 

In this section by applying the reconstructed action, we obtain the density contrast and 
corresponding power spectrum of the matter and compare it with that from the dark energy 



model. We start the procedure by reconstructing the screened mass function by substituting 
the numerical value for the first derivative of the action in Eq.( |3.18| ). The corresponding 
screened mass function for this action is plotted in Fig.(^). At redshifts above z > 1.5, 
the screened mass function approaches to unity in which the GR structure formation is 
recovered. It should be noted that the screened mass function does not depend only on 
the chosen MG action but also it depends on the wavenumber of the structures. One step 
further to indicate the probable different predictions of metric and Palatini formalism in 
matter density power spectrum we also plot the screened mass function in metric formalism 
Q"^^*{k,z) in Fig.(0) by applying the reconstructed dynamics of universe and MG action 
in Eq.(p8|). 




Figure 6: Screened mass function in Palatini formalism in terms of the redshift and wavenumbcrs 



of the structures for the action equivalent to the dark energy with the equation of state in Eq. (4.7) 



The comparison of two Figures @ and (0) shows the different scale dependence of 
screened mass function in both formalisms and consequently the different effects on the 
modification to the Poisson equations. A crucial point to indicate is that this difference arise 
from the structure formation effects. This is because the background dynamics H = H(z), 
is the same for two formalisms obtained independently from SNIa data. The difference 
comes from the definition of Ricci scalar R{g) = 6H + 12H^ and modified Friedmann 
Eq.( ^.39|) in metric formalism and also the screened mass function defined in Eq.( 3.3^ ). 



In Palatini formalism, in small wavenumbers, correspond to the larger structures, 
Q{k,z) descends for the later times, means weakening of the gravitational strength in the 
Poisson equation. On the other hand, in the smaller structures Q is larger than one (see. 
Fig.p) which makes a faster growth of the structures compared to the GR. The dependence 
of Q to the wavenumber is similar to the dispersion relation in optics in which different 
length structures' gravitational response is different The deviation of Q from unity 
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Figure 7: Screened mass function in metric formalism in terms of redshift and wavenumbers of 



the structures for the action equivalent to the dark energy with the equation of state in Eq. (4.7). 



and also its dependence on k affects the growth of the density contrast and subsequently 
the power spectrum of the structures. In Palatini formalism the deviation from Q = 1 is 
large at the larger scales, in contrast in metric formalism we always have a faster growth 
in structures and the the deviation is larger in small scales. These differences can put their 
fingerprint in LSS observations. In what follows we will solve the differential equation gov- 
erning the density contrast for various wavenumbers and finally obtain the corresponding 
power spectrum of the structures. 

In order to solve Eq.(3.29) we need to obtain the gravitational slip parameter. Sub- 
stituting the action in Eq.( p.25D , we reconstruct the gravitational slip parameter as shown 
in Fig.(|8|). The gravitational slip parameter is equal to one for higher redshifts and GR 
is recovered. For the later times at lower redshifts 7 deviates from GR. This deviation is 
larger for the small structures. The variation of Q and 7 with respect to GR indicates that 
the power spectrum of the structures should deviate in comparison with the DE models. 
It is worth to mention that the cosmic shear statistics as well as ISW effect is sensitive to 
the gravitational slip parameter and provides a promising tool for distinguishing between 
the models. 



For calculating the evolution of the density contrast we reexpress Eq.( 3.29| ) in terms 
of redshift as below: 



A(k, -h B(k, z)6k = 0, 



(4.13) 
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Figure 8: The gravitational slip parameter is plotted versus redshift for different wavenumbers of 
the structures. 



where the coefficients A and B are defined as 
A 
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'e ^ l + z 



l + z 



Qil + z) 9nl^ Q 



2k^ l + z 



'e 



.'E 



{l + z)Q — -Q -{l + z)Q )], 



(4.14) 



(4.15) 



and a{k,z) = 2A;V[2p + 9il^Q(l + z)], where k = k/Ho- For the case that Q' = g" ~ o 
and k > 9J1^Q(1 + z) we will recover the general form of density contrast evolution as: 



6 + 



E^ 

'e 



l + z 



(4.16) 



In this work within the Palatini formalism we solve the general form of Eq.( |4.13 ). The 
result for the reconstructed modified gravity is plotted in Fig.(|9|) for the various wave 
numbers. Here we take the initial condition of the density contrast from the Harrison- 
Zeldovich spectrum before the turn over point of the power spectrum around k ~ 0.02. 



By solving Eq.( 4.13|) we obtain the power spectrum of the structures at the present 



time as plotted it in Fig. (|T^ where the slope is nj = 1.21^019. It should be noted that at 
these large scales (i.e. k < 0.01) we don't have accurate data from the large scale structure 
observations ||5^. In order to compare the two models with nowadays LSS survey data, 
we should study the nonlinear regime of power spectrum which is not in the scope of this 
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Figure 9: The matter density contrast in Palatini formalism is plotted versus redshift for different 
wavenumbers. 
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Figure 10: The power spectrum of MG theory in terms of wavenumber(dashed line) is plotted in 
comparison to Harrison-Zeldovic power spectrum (solid line). For the MG model the slop of the 
power spectrum is 



1 91+0.13 
-■^•^^-o.ig- 



work. We note that the transformation of the linear power spectrum to nonlinear regime 
by analytical methods is strictly model dependent and it is not yet well understood p6[ . 

In order to compare power spectrum predicted from the equivalence DE and MG 
models, we obtain the density contrast in DE model for various wavenumbers at the present 



0.01 0.02 0.03 

k 



Figure 11: The fraction of MG to sDE power spectrum is plotted versus wavenumbers at the 
present epoch. 

time. We note that the only difference in the differential equation governing the density 
contrast evolution of the MG and DE is that, for the case of DE we let Q, 7 = 1. We plot the 
fraction of the power spectrum of the matter density in the MG to that in the equivalence 
DE in terms of wavenumber at present time in Fig.(^). While we have equivalent dynamics 
for the background expansion of the universe, the power spectrum from the MG and DE 
models are different. For the structures with the size of A; < 0.01 the MG power spectrum 
is smaller than the DE model, as the screened mass in these scales is less than unity. On 
the other hand for structures with the scale of k > 0.01 the screened mass is larger than 
unity and the ratio of the two power spectrums increases. The dependence of the power 
spectrum to the wavenumber can provide an observational tool to distinguish between these 
two models. Qualitatively we can also assert that as the screened mass function in small 
scales is bigger than unity, we will have an enhancement in mater density power spectrum 
in comparison to the equivalent sDE model. 

Another observational parameter in the structure formation is the dynamics of the 
growth index factor defined as / = din 6 /din a in terms of the redshift. We rewrite the 
differential Eq.( [4.13| ) in terms of growth index factor / as below: 

/ - + - Y^Jf - B{1 + z) = 0, (4.17) 



where A and B are defined from Eqs. ( [4.14| ) and ( |4.15 ). We plot the growth index in the 



MG and sDE in Fig.(^) with the observational data from the 2dF survey [^]. For the 
higher redshifts the growth index for the two models is the same while for the later times 
they diverge. An observational way to distinguish the MG from the DE model is measuring 
the growth index for various wavenumbers. For the case of DE model, the growth index 




Figure 12: The growth index function is plotted versus redshift for various wavenumbers. The 
observational data from the 2dF survey is represented with the corresponding error bar. 



is independent of the size of structure, as the structure formation equation for the scales 
larger than the Jeans length is independent of the k while in the MG model, screened 
mass relates the growth index of the structure to its size. Future surveys of the large scale 
structure may reveal the growth index in terms of wavenumber of the structures. 



5. ISW as probe to distinguish between DE and MG models 

In this section we examine the effect of MG on Integrated Sachs- Wolfe effect (ISW) [^8| 
and compare it with the dynamically equivalent dark energy model. ISW effect is one of 
the important cosmological probes which traces the effect of dynamics of the potentials on 
the photons of the cosmic microwave background radiation (CMB). The effect of MG on 



ISW effect has already been studied in several works in metric formalism as in |59|. This 
effect has been also studied in a class of /(i?)-gravities in Palatini formalism where the 
deviation from ACDM of this models occurs at a higher redshift (known as the early time 
f(R) gravities) ||6l|] and the results are compared with ordinary f{R) gravities. It was shown 
that the matter power spectrum is more sensitive to gravity modification than the ISW- 
effect, which is compatible to our results obtained in this section, although their specific 
choice of gravity action caused to a deviation in ISW effect in higher 1-moments [ |60[ . 

We start with the reconstruction of the perturbed potential $ which is a relevant 
parameter in the ISW and weak lensing observations. Using the Poisson equation of ( 3.16| ) 



for MG, $ is given in terms of the density contrast and screened mass function as follows: 
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Figure 13: The perturbed potential $ in MG tlreory in terms of redsliift for various wavenumbers. 
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Figure 14: The fraction of perturbed potential in MG to that in the sDE model in terms of redshift 
for various wavenumbers. 



Substituting the numerical values of matter density contrast and the screened mass 
function, we obtain the potential in terms of redshift for various wavenumbers. The result 
is plotted in Fig.(^). For comparison with the dark energy model, we express the relative 
magnitude of Potential in MG and sDE models as below: 



Q{k,z). (5.2) 



Fig.(|T^ shows this ratio in terms of redshift for various wavenumbers. Eq.(^[^) shows that 
this ratio is proportional to the relative magnitude of density contrast and screened mass 
function. In the zeroth order the density contrast growth is approximately equal to density 
contrast in sDE. The relative magnitude of potentials diverge from unity in large scales as 
the screened mass function k-dependance. Although the k-dependence of 7-parameter has 



a secondary effect on the growth of density contrast through Eq.( 4.13 ). 

Now we apply results from the dynamical evolution of potential in the ISW effect. ISW 
effect is caused by the time variation in the cosmic gravitational potential, as the CMB 
photons pass through the structures. In this effect the amount of the energy gain of the 
photons due to the gravitational potential of the structure may not be equal to the energy 
lose of the photons when they leave the potential well. The result of this effect is an extra 
anisotropy on the CMB map. 

In generic case, the CMB anisotropy due to the ISW is given by 



Qi{k,Vo] 



drje 



^{k,r]) -^{k,r]) ji{k{r]o-r])), 



(5.3) 



where 0; represents the 1-th moment of temperature contrast defined in Eq.(2.5), r is the 
optical depth of the photons which we can neglect it in ISW calculations, ji is the Bessel 
function and 77 is the conformal time and "." in this equation represents derivative with 
respect to the conformal time. Since in the smooth dark energy models <I> = — Eq. 
reduces to: 



e?^^(fc,r?o) 



d7] -2«>^^(fc,r/) ji{k{rjo - 7])) 



(5.4) 



superscript "sDE" on Q indicates ISW effect in sDE model. For calculating the ISW in 
the equivalent dark energy model, we use the numerical value for the $ which is obtained 
in the last section. On the other hand, since in the MG, the two potentials of ^ and ^ are 



not equal and relates by the slip parameter in Eq. (3.25), we eliminate in favor of <I> in 
Eq.(^.3|) and rewrite the ISW effect in the MG models as follows: 



qMG 
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z=0 
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dz 



(5.5) 



where the superscript of "MG" in 0; indicates the temperature anisotropy in the MG 
model. For simplicity in the numerical calculation, we reexpress time derivatives in terms 
of redshift and since $ is almost a constant function beyond the redshift of z ~ 1.7 (see 



Fig. 13), we put a cutoff in the numerical integration up to this redshift. 

In order to express our results in terms of observable parameter, we obtain the 0/ in 



terms of CMB power spectrum, Q. The variance of temperature fiuctuations, (5^^/ 6^ 
air 



( m 



I > can be related to 0; as below [29|: 
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(5.6) 



We note that for the small scale structures, the time derivative of the potential in 
Fig. (13) is almost zero. Hence we expect that in the ISW, small structures ( large Is ) 
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Figure 15: The difference of AC; — Cf^^ — C^'^^ between the equivalent DE and MG models in 
terms of I. 



will not be modified. So we examine only small moments for calculating deviation from 
the standard CMB power spectrum. Taking small Is in the range of / < 10 we calculate 
Cis both in the MG and the equivalent DE model. The numerical difference of the power 
spectrum ACi = C^^^ — Cf'^'^ between the DE and MG models in terms of / is shown in 
Fig.(|T5|). Here for the small Is we have more contrast than the large Is. 

Finally we express our analysis in CMB power spectrum for the MG and sDE models 
and compare them with the observational data. We use the CMB-fast code to plot the 
power spectrum for the equation of state defined in Eq.( [4.7D . As we showed in Fig. (p^), 
the difference in the power spectrum between the two equivalent DE and MG models is 
in the small Is and for the larger Is the power spectrum of the two models coincide to 
each other. For smaller Is, we use the numerical results of Fig. ( p!5| ) to obtain the power 
spectrum of the MG from the DE model. Fig.(^) compares the power spectrum from the 
MG and DE model with the latest WMAP-7 years data [32|. We note that this deviation 
is smaller than the cosmic variance in the WMAP data. 



6. Conclusion 

One of the important questions in the problem of the acceleration of the universe is either 
this acceleration is produced by a dark energy fluid or that is a manifestation of the 
gravity modification. It has been shown that the expansion history of the universe can not 
dynamically distinguish between these two models and at this stage they are equivalent. 

In this work we compared various aspects of the structure formation in the Modified 
Gravity (MG) and smooth Dark Energy (sDE) models. We used an ansatz for the equation 
of state of a dark energy model and obtained the equivalent modified gravity . Using a 
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Figure 16: Power Spectrum of CMB for Modified gravity and equivalent dark energy model in the 
Palatini formalism. 

synthetic SNIa data set from the SNAP observation, we did the inverse problem approach 
and obtained the effective action of the gravity (i.e. f{R)). It has been shown that the 
structure formation unhke to the background dynamics can distinguish between these two 
models. 

In the modified gravity, the modification of the Einstein-Boltzman equation for the 
structure formation has extra terms compare to the sDE model. This modification can be 
given in terms of screened mass function and gravitational slip parameter. We obtained 
these two parameters in the MG model and solved the differential equation for the evolution 
of the structure formation. Using the Harrison-Zeldovich initial condition for the structures, 
finally we obtained the power spectrum and the growth index of the structures in the two 
scenarios in the linear regime of the structure formation. We have shown that for the 
structures with the larger size, the effective spectral index at the present time is slightly 
more than one. Prom the observational point of view, sampling of large scale structure more 
than lOOMpc is needed to test this effect. On the other hand, we showed that the growth 
index parameter in the MG unlike to the sDE models is a scale dependent parameter. 
Next generation of peculiar velocity survey may measure this effect. Pinally we obtained 
the effect of the MG structure formation on the Integrated Sachs- Wolfe effect on CMB 
map. We showed that the deviation of the CMB power spectrum due to the MG is larger 
in the small Is, however it is not too large to distinguish between the two models due to 
the uncertainty from the cosmic variance. 

A. Structure formation probes in metric formalism 




Here in this appendix, we obtain the evolution of matter density contrast in the metric 



formalism and consequently derive the matter power spectrum and growth index in order 
to compare them with that in the Palatini formalism. 
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Figure 17: The gravitational slip parameter in metric formalism for different wavenumbers. The 
Long dashed line correspond to small scales k — 0.02, where we have the most deviation in this 
parameter from GR value. 



In order to determine the matter density evolution we use Eq.(3.37), where we can 
rewrite it, in terms of redshift as below: 
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(A.l) 



The above equation has the same form as in the Palatini formalism, Eq.( [4.16 ), however the 
difference is in the value of the screened mass function and the gravitation slip parameter 
in the two modified gravity formalism. 

The screened mass function is obtained in Sec.(Q) and plotted in Fig.(^. For calculat- 
ing the gravitational slip parameter, due to the slow change of the action over the Hubble 
time we neglect the term of 5'dF{R) in the perturbation of the trace of the field equation 



(13.331) . We follow the same procedure an in the Palatini formalism and obtain the gravi- 
tational slip parameter in terms of reconstructed F from Eq. ( |3.4C ) and the screened mass 
Qmet from Eq.( |3.38|) which results in 
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Figure 18: The density contrast of the mater in metric formalism in terms of redshift and wavenum- 
ber is plotted, parameter in metric formalism for different wavenumbers. 



The gravitational slip parameter in the metric formalism is Plotted in Fig.(|l7|) for differ- 
ent wavenumbers. Similar to the Palatini case, Fig.(^, the gravitational slip parameter 
deviates from its GR value 7 = 1 in small scales. Knowing the screened mass function and 
gravitational slip parameter, we solve the differential Eq.( [A.f| ) governing the evolution of 
density contrast. The results is plotted in Fig.(p^). 

In the next step we compare the relative magnitude of power spectrum of metric-MG 
to sDE, Puet/PsDE and Palatini-MG to sDE , Ppai/PsDE both plotted in Fig.®. As we 
except because the variation of screened mass function from GR value Q > 1 in small scales, 
we get a amplification in power spectrum in large A:'s, which shows a different prediction 
than the Palatini formalism in which there is a suppression in large scales. 

Using Eq.( [A.l| ) we can also reexpress the growth index function / = ^}"^ in terms of 
redshift as: 

f + l^-^l^ltQmet{l + Zf = ^ (A.3) 

For comparison with the Palatini formalism, we plot the growth index for large and small 
scales in both metric and Palatini formalism in Fig.(^). 
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Finally the last comparison is done for the relative magnitude of gravitational potential 
^ in both formalism. From Poison equation, it is straightforward to calculate this ratio as: 

"'"'■fe-j^S^ (A.4) 



Met 



1.4 

1.2 
1 

0.8 
0.6 
0.4 
0.2 

0.005 , 0.01 0.015 0.02 

k (h Mpc) ' 

Figure 19: The relative power Spectrum of Modified gravity and smooth Dark Energy is plotted 
versus wavenumber for Metric formalism(solid line), Palatini Formalism (dash-dot line) 

The results for three wavenumbers is plotted in Fig.(|2l]). The relative magnitude of gravita- 
tional potential in higher-redshifts converge to unity for all wavenumbers. This is because 
the screened mass functions, in both formalism converge to unity and the evolution of 
matter density contrast in deep matter dominated era is proportional to the scale factor. 
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